A New Approach to the Calculation of Some properties of 

Vacuum and 7r, a Mesons 



Hong-shi Zong 1,2 , Xiao-hua Wu 1,3 , Xiao-fu Lii 3 , Chao-hsi Chang 2,4 and En-guang Zhao 2,4 

1 Department of Physics, Nanjing University, Nanjing 210093, P. R. China 

2 CCASTfWorld Laboratory), P.O. Box 8730, Beijing 100080, P. R. China 

3 Department of Physics, Sichuan University, Chengdu 610064, P- R- China 



(N 
Or 

m 



4 Institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing 100080, P. R. China 



Abstract 



> 
(N . 

It is shown in a systematic way that the instanton dilute liquid model can 
■ be combined with the global color symmetry model to form a new nonper- 

o ' 

turbative QCD model. Based on the quark propagator derived in the instan- 
r~| ■ ton dilute liquid approximation, the quark condensate (qq), the mixed quark 

O ■ gluon condensate gsiqG^a^q) , the four quark condensate (qTqqTq) and ten- 

sor, pion vacuum susceptibilities have been calculated at the mean field level 
V"! ■ in the framework of the new nonperturbative QCD model. The numerical 

results are compatible with the values obtained within other nonperturbative 
approaches. The calculated masses and decay constants of ir and a mesons 
are also close to the experimental values. 
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I. INTRODUCTION 



The non-perturbative structure of the QCD vacuum are characterized by the vari- 
ous condensates: such as the quark condensate (qq), the mixed quark gluon condensate 
gsiqG^a^q), the four quark condensate (qTqqTq), etc. These condensates would vanish in 
a perturbative vacuum but do not vanish in the non-perturbative QCD vacuum, and are 
of essential for describing the physics of strong interaction[l,2]. Meanwhile susceptibilities 
of vacuum are also important quantities of strong interaction physics. They directly enter 
in the determination of hadron properties in the QCD sum rule external field approach[3- 
5]. In particular, tensor susceptibilities of the vacuum[6] are relevant for the determination 
of the tensor charge of the nucleon which is connected through deep-inelastic sum rules 
to the leading-twist nucleon transversity distribution [7]. The strong and parity-violating 
pion-nucleon coupling depends crucially upon x n , the n susceptibility [8]. In addition, as 
Goldstone bosons, the n and a mesons provide ideal samples of exploring models of hadron 
structure in terms of the elementary degrees of freedom in quantum chromodynamics(QCD). 
They are then the ideal objects in investigating the nonperturbative aspect of QCD. Up 
to now, the main techniques for studying the nonperturbative aspects of QCD are lattice 
gauge theory [9], QCD sum rules[l-5], the instanton approximation[10-12] and some other 
phenomenological approaches. They have achieved quite success in describing different char- 
acteristics of low-energy QCD. However, each of these approach has its strengths and weak- 
nesses. It is interesting to combine two kinds of different nonperturbative QCD approaches 
to determine the vacuum condensates , vacuum susceptibilities and some properties of 7r 
and a mesons within a unified model. 

In the present paper, we try to address the above issues by combining the instanton dilute 
liquid model and the global color symmetry model(GCM) in a systematic way. Recently 
there are evidences that GCM provides a successful description of various nonperturbative 
aspect of strong interaction physics and the QCD vacuum as well as hadronic phenomena 
at low energies[13-15]. 
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Although the instanton vacuum has been shown to provide a good phenomenological 
description of many hadronic properties [11, 12], it has also some weaknesses[16]. To be con- 
crete, it assert the vector self energy function of quark properator in instanton medium 
A(q 2 )=l which is not consistent with the corresponding one in Lattice QCD[17]. In addi- 
tion, the q 2 3> lGeV 2 behavior of scalar self energy function B(q 2 ) of quark propagator 
in instanton medium disagrees with QCD. Therefore, how well does the instanton dilute 
liquid vacuum model the QCD vacuum is an interesting question. Furthermore, it is also 
an interesting question to what extent can the inclusion of perturbative one gluon exchange 
term improve the predictions of the instanton model in a systematic way. 

To proceed with our model one just explores its consequences, applying it to different 
systems , where possible, by comparisons with experimental data and other nonperturbative 
QCD approaches on their common domain. In this way we try to understand how well does 
our nonperturbative QCD model work. In order to keep this paper self-contained and easy 
to follow, a brief introduction of this nonperturbative QCD model is described in Sec. II. 
Some properties of the ir and a mesons, such as f n , m f , f a and m a , will also be evaluated 
in Sec. II. From the determined generating functional, the vacuum condensates and vacuum 
susceptibilities are calculated in Sec. III. Finally a brief discussion is given in Sec. IV. 

II. A NONPERTURBATIVE QCD MODEL AND SOME 
PROPERTIES OF tt and a 

It is well known that the QCD vacuum has topological charge because QCD is a non- 
abelian gauge theory. In order to consider the contribution of the nonperturbative QCD 
vacuum , instantons have to be introduced. It is believed that the instanton dilute liquid 
vacuum is a good approximation to the QCD vacuum configuration[18]. The perturbative 
contribution is deduced from the expansion about the instanton dilute liquid configuration, 
i.e., = Aft + A^, here A^ are the fields of the instanton liquid and A ^ are the quantum 
fluctuation fields of gluons. We will show in a reasonable way that the effect of the instanton 
liquid can be replaced by the Goldstone bosons. At the same time the hard perturbative 
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term, in principle, can be calculated in this way. Thus there is no double counting in our 
approach. 

Now we will use the path integral method to consider the contribution of the instanton 
dilute liquid. The generating functional for QCD in the Euclidean metric is given by 



Z[J,fj,ri] = J VqDqVAexp{-S + J d 4 x(fjq + q V + J^)J , 



(1) 



where 



n I r<a fta I 

H T (IV^ (JLV | i 



S = J d A x^q 

and Gp V = d^A a v — d u A a tl + gf abc A b ^Al. We leave the gauge fixing term, the ghost field 



term and its integration measure to be understood. It is believed that the dilute liquid of 
instantons is a good approximation to the extreme of the pure gluon fields action[18]. The 
pure gluon fields action may be expanded about the dilute liquid of instantons in terms of 
the quantum fluctuation variables A°^ . With the functional expansion, the pure gluon fields 
action can be written as follows 



x 



= I ^x\g%{x)G^{x) + \ I d*xd*yA^(x)Vf u (x,y)A b J>(y) + 



^ (tv \ / (iv v / 1 2 J & n \ / (IV \ T o / - -v \& / 1 ' 

(2) 

here Gf u = d^A a J - d v A* + gf^A^A^ 1 , and the term V ab v (x, y) can be regarded as the two 
point function of the quantum fluctuation fields. 

Substituting Eq.(2) into Eq.(l), the generating functional for QCD reads 



Z[J,fj, V } = JVqDqVAiexpl-S l + J d 4 x i gq^A^q + q V + r ] q + J a ^ + j;A a ; 



X -J d 4 xd 4 yA^(x)V^(x,y)A b J(y) 



+ 



(3) 



where 



S l = Jd'x^q ^[d^-ig^A^ 



4 - (IV" (IV ■ 

Using functional differentiation technique, Eq.(3) can be written as follows 



4 



Z[J,ri,rj\ 

= J VA«exp {-/ d 4 x(-^ig^A^^ - J^) -\f d'xd 4 yA^(x)Vf u (x : y)A^(y) + ■ 
x f VqVqexp {-^ + | d 4 x(w + q V + J^)} ■ 

Eq.(4) can be used to separates out the contribution of the instanton dilute liquid from 
the quantum fluctuation field. Thus we can study the effects of instanton dilute liquid and 
its quantum fluctuation fields separately. Now, we will pay attention the contribution of 
the instanton dilute liquid. In addition, it should be noted that it is interesting to consider 
further the effect of the perturbative gluon two-point function to the prediction of instanton 
model(or vice versa). 

Using again the functional differentiation technique, the generating functional for the 
contribution of the dilute liquid of instantons can be written as follows 

Z l [J,fj,v] = f VqVqexp {-S l + f d 4 x(fjq + qq + J^A?)} 

/if X a d 1 

VqVqexp j J d 4 x(-q-f^d^q + fjq + qr] + igq— 7/^7^) j exp (W(J)) , (5) 

where 

w[J] = f d'x {--g^ + . 

It is easy to see that Eq.(5) describes the interaction between the quark and the dilute liquid 
of instanton. 

It is obvious that exp(W(J)) can be expanded in powers of J® as follows 

ex7mJ=\) = I ^vl^'VKWM + W R [J] (6) 
where W R [Jfi = E^°= 3 h. I dx ^ ' ' ' dx n D X' a »n( Xl ' " ' ' x n) U i=i J ^( x i) , and according to the 
model of the instanton dilute liquid, D^ v (x,y) = (0 | A"j(x)A^(y) | 0) is evaluated on the 
vacuum configuration of the instanton dilute liquid. D°^ v is directly related to the gluon 
condensates. 

Since exp(W[J = 0]) is a normalization factor, it can be omitted for convenience in 
the following. Introducing an auxiliary bilocal field B e (x,y) as in [13-15], the generating 
functional of Z l [J,r),fj] can be written as 



Z l [0, 77, 77] = exp 



X 



£77(2;) 2 Sr](x) 

J VqVqVB e (x,y)ex V {-S{q,q,B e (x,y)] + J d 4 x(fjq + qv)} , (7) 



and 



S[q,q,B e (x,y)] = J J d A xd 4 y 



, r 0n , v B 9 (x,y)B e (y,x) 

q(x)G (x,y-[B])q(y)+ 2gW{x _ y) 



with 



G" 1 ^, y; = 7 • - y) + V# V, y) , 



(8) 



where A e = K a C b F c is determined by Fierz transformation in color, flavor and Lorentz space. 
Since we have to take the average in the color space for calculation of the g 2 D ab (x — y), 



we have choosed a suitable gauge B ab v [x — y) = 5^5 D(x — y) in deriving the action 
S[q,q,B 9 (x,y)\. 

Neglecting W0?[J^] as in [13-15] and performing the functional integral over Vq and Vq 
in Eq. (7), we obtain the modified GCM generating functional as 



Z G cm[v, v] = J T>B e (x, y) exp(-%, 77, B 9 (x, y)]) , 



(9) 



where 



S[fj,ri,B (x,y)] = -Trln 
+ 



p6(x-y) + -A e B e (x,y) 



B 9 (x,y)B 9 (y,x) 



+ fj(x)G(x,y; [B e }) v (y) 



(10) 



2g 2 D(x-y) 

From Eq.(10), it is easy to see that the S[fj, 77, B e (x, y)\ is the action of the bilocal bosonized 
fields. Here bilocal fields arise from the interaction between the quark fields and the dilute 
liquid of instanton. 



The saddle-point of the action is defined as 5S[T], 77, B 9 (x, y)]/5B 9 (x, y) 
given by 



»7=»7=0 



= and is 



B e (x -y) = g 2 D{x - y)tr^c[A e G (x - y)], 



(11) 



where G stands for G[Bq]. These configurations are related to vacuum condensates and 
provide self-energy dressing of the quarks through the definition: 

E(p) = l^B^p) = J Jxe** [lA 9 B e (x)\ = * 7 • p[A(p 2 ) - 1] + B(p 2 ), (12) 

where A(p 2 ) and B(p 2 ) are the scalar and the vector self energy function of quark propagator, 
respectively. This dressing comprises the notion of "constituent quarks" by providing a mass 
M{p 2 ) = B(p 2 )/A(p 2 ), reflecting a vacuum configuration with dynamically broken chiral 
symmetry. 

Recalling Eq.(8), at the mean field level(ie., the bilocal field B®(x,y) is simple replaced 
by the bilocal field vacuum configuration B^(x — y) in Eq.(8), physically this means that we 
only consider the influence of instanton medium to the quark propagator), we get 

(0\T[q(x)q(y)]\0)- 1 = 7 • dS^(x - y) + -jB°(x - y), (13) 

with this equation the vacuum configuration of the bilocal field Bq{x — y) is related to 
the quark propagator. It is well known that the propagation of light quarks was formulated 
quantitatively by treating quarks propagating in the instanton medium[12] and refined in its 
treatment of nonzero-mode propagation by Pobylitsa[19]. In momentum space, The quark 
propagator derived in the instanton dilute liquid model can be written as[19]: 



° W q 2 + m 2 (q 2 ) W > V ™ c ^ J(2n)~Wkk 2 ^(k) ^ ' 

where <p(q) = 7rR 2 £[I (z)K (z) - I^K^z)], in which z = M5. I n (z) (K n (z)) 
(n = 0, 1) are the first (second) kind modified Bessel functions of order n. p ~ (200 MeV) 4 
is the average density of the instantons. R = | fm is the average radius of the instantons. 

As similar to the local field, the quantum bilocal fields B q (x,y) can also be expanded at 
the saddle point B®(x — y) and the fluctuations with respect to B®(x — y). Each fluctuation 
can be separated into two factors: the motion of the center of mass and the relative motion. 
A separation of the internal and center-of-mass dynamics is achieved by considering the 
normal modes of the free kinetic operator of the bilocal fields in a manner that is analogous 



to the interaction representation of standard quantum field theory. As only the excitation 
of the Goldstone bosons is taken into account, we showed that the vacuum configuration of 
the bilocal field B(q 2 ) and the internal motion of the Goldstone bosons are equivalent to 
the dynamical mass m(q 2 ) of quarks in the instanton dilute liquid approximation [20]. The 
mesons 7r and a generated by the bilocal field in our model is then determined as 



where the 6 takes only n, a and B{x — y) is the Fourier transformation of the m(q 2 ). 

So far, we have finished the work of combining the instanton dilute liquid model and 
the GCM in a systematic way. As every parameter has been fixed in the instanton dilute 
liquid approximation, there are no any free parameter in our nonperturbative QCD model. 
It should be noted that strong interaction phenomena provide evidences in favor of the 
instanton structure of the QCD vacuum. Lattice QCD [18,21] calculation shows that the 
instanton dilute liquid is a good approximation of QCD vacuum. In addition, we want to 
stress that the vacuum configuration of the bilocal field in GCM and the internal motion of 
the Goldstone bosons are directly related to the instanton dilute liquid approximation of the 
QCD vacuum. The results given below then reflected directly how well does the instanton 
dilute liquid vacuum model the QCD vacuum. 

Based on Eqs.(10) and (15), the effective action of the Goldstone particles can be written 

as 



With Eq.(16), following the same way as shown in Refs. [13,20], we obtain the normal- 
ization constants and the masses of 7r and a mesons as 



B°(x, y) = B e (x - y) + E 4> 9 q (-^ JL )B(x - y), 



(15) 



e 




(16) 



J TT 




OO 



sds 



W 2 + 2B s<w + gB 3^B + 3sB 2 (^) 2 sB 2 (l + 25f f 



(17) 



(s + B 2 ) 2 (s + 5 2 ) 3 
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9 3m f 00 , B(s) 

m *=2^/„ s<fa 7Tl4)' (18) 



/ 2 = — f 

8jr' 2 Jo 



OO 

sds 



l + 25f + ^-«(f) ; 



(s + 5 2 ) 



2 



5 2 
+ 



2 



1 + 25f + s50 + s(f ) 2 j 5 2 ( S - B 2 ) [l + 25(^?) 2 



ds 



sB 2 (s-B 2 



(s + B 2 ) 3 (s + B 2 Y 



(19) 



3 f°° B A (s) 

ml = W;i sds ¥TWW' (20) 

With Eqs.(14), (17), (19) and (20), we obtain the normalization constants 7^=86 MeV, 
/ cr =77 MeV and m cr =649 MeV. As the n normalization constant can be simply related to the 
decay constant of n when the vector self-energy function of quark propagator A(q 2 ) = 1[22], 
the calculated result of f w is quite close to the experimental data ^=93 MeV. Meanwhile 
the obtained m a is quite reasonable which is in the scope of experiments(500-700) MeV[23]. 
Furthermore, the expression of m n reproduces a PCAC result that ^ only if the 
current quark mass m/0. This agrees well with the Goldstone-Boson feature of the pion. 
Moreover, we fit m T = 139 MeV as the current quark mass m is taken to be 3.7 MeV. 

From Eq.(18) and the calculation of the two quark condensate (qq) below, we have 

fl m l = -2m(gg). 

which is the current algebra result derived by Gell-Mann, Oakes and Renner(GMOR). 

III. CALCULATION OF SOME PROPERTIES OF VACUUM 
A. EVALUATION OF VACUUM CONDENSATES 

Once the generating functional in our model are determined, one can calculate the (qq), 
gsiqG^a^q) and (qTqqTq) vacuum condensate at the mean field level. Previous studies of 
these condensates include QCD sum rules where it was treated as fitted parameter in the 
analysis of various hadron channels [1,2], quenched lattice QCD [24], the instanton liquid 



model [25] and the model of confining gluon propagator used in the Schwinger-Dyson form 
for GCM[26]. 

With the modified GCM generating functional, it is straightforward to calculate the 
vacuum expectation value of any quark operator of the form 

O n = {q^%q a ){q j2 K%q a ) ■ ■ ■ (q jn A%\ n q in ), (21) 

in the mean field vacuum. Here the AW stands for an operator in Dirac, flavor and color 
space. 

Take the appropriate number of derivatives with respect to external sources rji and rjj of 
Eq.(9) and putting rji = fjj = 0[27], one obtains 

(: O n :) = (-)" Xl[( _ ) PA iHl ' ' ' A i"L(G'o)ilip(l) ' • • (G ) in jp(n)], (22) 

p 

where p stands for a permutation of the n indices. In particular we obtain the nonlocal 
quark condensate (: q(x)q(0) :) as 

(: q(x)q(0) :) u = (-)tr, c [G (x,0)} 

= ( - AN ) r ^3 e w 

1 c) Jo {27cy P 2 A 2 {p 2 ) + B 2 {p 2 ) 

- (-)^- ["dec B{S) [2 Jl ^ l (23) 
" ( } 16n 2 Jo ^ sA%s) + B 2 (s) [2 ^ J' W 

where \i is the cutoff which we choose to be 1 GeV 2 as in Refs. [26] and [28] . It is natural to 
introduce a cutoff because our model is an effective theory and one expects that there will 
be a characteristic scale in which the theory is effective. This situation is very similar to the 
determination of vacuum condensates in Nambu and Jona-Lasinio(NJL) model. 
At x=0 the expression for the local condensate (: qq :) is recovered: 

A :>, = H*r,c[G„(*,0)]U = H ±, £ dss -^L^ . (24) 

The nonlocality g(x 2 ) can be obtained immediately by dividing Eqs.(23)by Eq.(24). In 
Table I, we display the nonlocal quark condensate g(x 2 ) in the text. By comparing it with 
the corresponding results mentioned in Ref.[29], we demonstrate that the nonlocal quark 
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condensate is very robust with respect to using our approach or the models of confining 
gluon propagators. 

Table. I. The nonlocal quark condensate g(x) = (: q(x)q(0) :)/(: q(0)q(0) :). 
x 2 (GeV- 2 ) 0.0 2.0 4.0 6.0 8.0 10 12 14 16 18 20 
g(x 2 ) 1.0 0.89 0.79 0.71 0.63 0.55 0.50 0.44 0.38 0.34 0.30 

Another important vacuum condensate is the nonlocal four quark condensate in the mean 
filed vacuum. From Eq.(22), one has 

(: q(x)A^q(x)q(y)A^q(y) % 
= -tr jC [Go(y,x)A^G (x,y)A^} + tr lC [G (x,x)A^}tr 7 c[G (y,y)A^}. (25) 

By means of Eq.(25), one can calculate all kinds of nonlocal four quark condensates at 
the mean field level. For instance in case of A* 1 ) =A^='y fl ^- one find from Eq.(25): 

9(^)7^9(^)5(0)7^9(0) 0,. 



o Jo (27r) 4 (27r) 4 

Mf) A(q 2 ) 

A 2 (p 2 )p 2 + B 2 {p 2 ) A 2 {q 2 )q 2 + B 2 {q 2 ) 



4 3 B(p 2 ) B(q 2 ) 



d 4 p d 4 q 



A 2 (p 2 )p 2 + B 2 (p 2 ) A 2 (q 2 )q 2 + B 2 (q 2 ) 



(26) 



Similarly at x=0 the expression for the local four quark condensate (: q'Jfj.-^-qq'J^-^-q '■) is 



recovered: 



v / ^xrf d 4 p B(p 2 ) l2 . ,4. . 2 

<= tr,f *r,f 9 % = M 3 )[/ o ^ W % V) ] 2 = H 5 <= «9 2 , (27) 

i.e. for the local four quark condensate, our result is consistent with the vacuum saturation 
assumption of Ref. [1]. However, if one considers the nonlocal four quark condensate, it 
should be noted that the contribution of the second term of the right-hand of Eq.(26) is 
about a 30% correction and can not be neglected [30]. 

The introduction of nonlocal condensates [31] have been shown to be useful for repre- 
senting the bilocal vacuum matrix elements needed for the pion wave function [32] and pion 
form factor [33] for low and medium momentum transfer. In the external field method of 
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QCD sum rules, one does not carry out OPE for the vacuum matrix element of the bilocal 
operators but introduce new phenomenological parameters to characterize the space-time 
structure of the nonlocal condensates. Previous studies of calculating this nonlocal conden- 
sate include the instanton liquid model [34] and the model of confining gluon propagator used 
in the Schwinger-Dyson form for GCM[29]. Although (: q(x)q(0) :) is not gauge invariant, 
it is useful for describing the scalar part of nonperturbative quark propagator. Therefore, 
we use our modified version of GCM to calculate the (: q(x)q(0) :) and nonlocal four quark 
condensate. 

As far as the mixed condensate g s (qG 'q) , one can use the method described by 
Refs.[26] to obtain the mixed condensate in Euclidean space 

/ „ U v k / u^u 27 B[2A(A - l)s + B 2 ] io rn , 2 B(2-A). /ooN 

9 .(5G^>„=(-) (l ^){ T / o dss L_l_ l + 12 ^ ^J^} (28) 

In Table II. we display the results for (qq) and g s (qG 'q) in our approach and compare 

it with the corresponding values which were obtained by other nonperturbative approaches: 

QCD sum rules [3], quenched lattice QCD [24], the instanton liquid model [25] and the model 

of confining gluon propagator used in the Schwinger-Dyson form for GCM [26]. 

Table. II. (qq) and (g s qG q) in different non-perturbative approaches 

{-)(qq)s (-)(gsqvGq)* 

this work 217 MeV 429 MeV 

QCD sum rules[3] 210 -230 MeV 375 -395 MeV 

quenched lattice [24] 225 MeV 402 -429 MeV 

instanton liquid model [25] 272 MeV 490 MeV 

the model of confining gluon propagator in GCM [26] 150 -180 MeV 400 -460 MeV 

Table II. shows that our results for (qq) and gs^G^a^q) are compatible with the 
values obtained within other nonperturbative methods, especially in QCD sum rules[3] and 
quenched lattice QCD [24]. 

B. EVALUATION OF VACUUM SUSCEPTIBILITIES 
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Now we pass to the calculation of pion and tensor vacuum susceptibilities. Quite different 
values of tensor vacuum susceptibility have been reported in literature [35-38]. Therefore it 
is interesting to address this issue in our model. 

In the external field of QCD sum rule two-point method, one often encounters the quark 
propagator in the presence of the J r (y) = q(y)Tq(y) current (r stands for the appropriate 
combination of Dirac and flavor matrices). 

S c £{x) = (0\T[q c a (x)q^o)}0)jr 

= S%^ T {x) + S%^ NP {x) ) (29) 

where Sg ,PT (x) is the quark propagator coupled perturbatively to the current and Sg ,NP (x) 
is the nonperturbative quark propagator in the presence of the external current J r (one 
should note that the external current should be taken to be J r <pr, where 0r is the value of 
the external field. In what follows, to simplify the notation, we will take the 0r=l, which 
does not affect our results). The vacuum susceptibility x T i n the QCD sum rule two-point 
external field treatment can be defined as [39] 

S$ r * p (x) = {0\:&(x)<${0):\0)^ 

= —^T a p5 cc , X r H(x)(0\ : g(0)g(0) : |0>, (30) 

where the phenomenological function H(x) represents the nonlocality of the two quark non- 
local condensate. Note that H(0) = 1. 

The presence of external field implies that S^ r (x) is evaluated with an additional term 
AL = —J T ■ 0r added to the usual QCD Lagrangian. In three point method of QCD sum 
rule, if one takes only a linear external field approximation, the S^ /3 r ' NP (x) in Euclidean 
space is given by 

S:f' NP (x) = J d*y e-^(0| : g^)g%)r g %)g-'(0) : |0>. (31) 
Using Eqs.(30), Eq.(31) and note that H(0) = 1 we have 
- ^S cc ,T^ X r (0\ : 5(0)9(0) : |0> = J d% e"^(0| : q c a (0)q e (y)m y )q^0) : |0>. (32) 
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Multiplying Eq.(32) by Tp a 5 cd , we get 
16tt 2 



r 

X a = 



tr 7 (rr) / ^ e ~ tq ' y{ ° 1 ■ ^)Tq c m\y)T q \y) : |0), (33) 
with a=-(27r) 2 (0| : q(0)q(0) : |0). Eq.(33) shows that the vacuum susceptibilities originates 
from the nonlocal four quark condensate contribution. This conclusion is the same as that 
of Ref.[38] which addressed the problem from a completely different viewpoint, using the 
concept of duality (more details can be found in Ref.[38]). 

In the case of tensor current (r = a^), from Eqs.(26), Eq.(33) and the fact tr^a^a^) = 
48, we have the tensor vacuum susceptibility x Ta: 

T 

X a 

= Lrm q ^)jdye '^[/p^^e — — - e 



rtJ- 

Lim„^ 3 / sds 
Jo 



B(s) 



2 

2 



0.23 GeV% 



_A 2 s + B 2 (s)_ 

with x T(1 an opposite sign and a factor (4 x tt 2 ) larger than the definition in Refs. [35,36,38]. 
In order to compare it with the estimation in Refs. [35,36,38], we have 

T 

X " - - 0.0058 GeV 2 , (35) 



(-4 X 7T 2 ) 

which is very close to the result -0.0055 GeV 2 obtained in Ref.[38]. Notice, that the result 
of present paper is less than the estimations -0.010 ~ -0.008 GeV 2 obtained in Ref.[36], and 
has opposite sign and is several times larger than the previous estimation in Ref. [35]. In 
order to make it easy to compare our result with previous estimations of vacuum tensor 
susceptibilities, a simple table is listed below 

Table. III. Vacuum tensor susceptibilities(Gel /2 ). 



Ref[35] Ref[36] Ref[37] Ref[38] this work 
+0.002 -0.008 +0.009^+0.017-0.0055 -0.0058 

Now we turn to the calculation of Pion vacuum susceptibility(The Pion vacuum suscep- 
tibility is crucial to the strong and party- violating pion-nucleon coupling). In the case of 
pseudoscalar current, from Eqs.(26) and (33), we have the Pion vacuum susceptibility x^a: 
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- i^r y d y tr 7 |y ^ AV + jB2(p2) e Tsj^ ^2 + 52(^2) e 7s 
= 3 r sds ~r, Kwrr = 2-56 GeV 2 , (36) 

which is compatible with the range y^a ~ (1.7 — 3.0) GeV 2 obtained within a phenomeno- 
logical approach[39]. The integral in Eq.(36) is divergent, as is inherent in effective models, 
hence a cutoff is introduced in the integral of Eq.(33) as in the Nambu and Jona-Lasinio 
model. 

IV. Discussion 

In the present paper, we first gave a systematic way to combine two kinds of different 
nonperturbative QCD approaches, ie., the instanton dilute liquid model and the GCM, then 
based on the quark propagator derived in the instanton dilute liquid approximation, we have 
determined the quark condensate (qq), the mixed quark gluon condensate gsiqG^a^q) , the 
four quark condensate (qTqqTq) and tensor, pion vacuum susceptibilities at the mean field 
level in the framework of our nonperturbative QCD model which used the quark propagator 
derived in instanton dilute liquid vacuum as input. The numerical calculation shows our 
results is compatible with the range obtained within other nonperturbative approaches. In 
particular, the numerical calculation of the tensor vacuum susceptibility shows that our 
results is very close to the one obtained in Ref . [38] and has opposite sign to that of Refs. [35] 
and [37]. As has been pointed out in Ref. [35] experimental measurement of the tensor charge 
of the nucleon is possible, and one might be able thereby to test the theoretical predictions 
of the tensor susceptibility in the future. Meanwhile we have also calculated the masses and 
decay constants of n and a mesons. The calculated results agree with experimental data 

It should be noted that the cutoff must be introduced in the calculation of the vacuum 
condensates and susceptibilities. This is consistent with the philosophy of the model: it is 
applicable below the chiral symmetry breaking scale(about 1 GeV). Other nonperturbative 
QCD models, such as the NJL model and chiral quark model has the same limitation. Some 
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quantities, such as the mass and normalization constants of ix and a, are not sensitive to 
this cutoff. For example, if one take the same cutoff /x = lGeV 2 in the integral of Eq.(17) 
and Eq.(19), one get f n = 85 MeV and f a = 76 MeV which is scarcely affected by the cutoff. 
This is because every quantum fluctuation in GCM can be separated into the motion of the 
center of mass and the relative motion. The relative motion factor(form factor) plays the 
role of suppressing the large momentum behavior in the integral. However if one calculate 
the vacuum properties(such as vacuum condensates and susceptibilities), there are no such 
form factor to suppress the large momentum in the integral(some integral such as Eq.(36) 
is divergent), hence a cutoff is necessary in calculating the vacuum properties. 

Finally, we like to emphasize that there are no any free parameter in our model, the only 
input is the quark propagator in instanton vacuum. With this model, the calculated results 
of masses and the decay constants of the mesons n and a and some vacuum properties agree 
with experimental or experiential dates quite well. We conclude that quark propagator in 
instanton medium contain valuable information of nonperturbative QCD and is consistent 
with recent models of confining gluon propagators used in the Schwinger-Dyson formalism 
for the global color symmetry model. 
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